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Abstract 

We have derived long series expansions of the percolation probability for site 
and bond percolation on directed square and honeycomb lattices. For the square 
bond problem we have extended the series from 41 terms to 54, for the square 
site problem from 16 terms to 37, and for the honeycomb bond problem from 
13 terms to 36. Analysis of the series clearly shows that the critical exponent 
(3 is the same for all the problems confirming expectations of universality. For 
the critical probability and exponent we find in the square bond case, qc = 
0.3552994 ± 0.0000010, /3 = 0.27643 ± 0.00010, in the square site case qc = 
0.294515 ± 0.000005, (3 = 0.2763 ± 0.0003, and in the honeycomb bond case 
qc = 0.177143 ± 0.000002, j3 = 0.2763 ± 0.0002. In addition we have obtained 
accurate estimates for the critical amplitudes. In all cases we find that the 
leading correction to scaling term is analytic, i.e., the confluent exponent A = 1. 



PACS numbers: 05.50.+q, 02.50.-r, 05.70.Ln 



1 Introduction 



Directed percolation (DP) was originally introduced by Broadbent and Hammersley 
(1957) as a model of fluid-flow through a random medium and has since been associ- 
ated with a wide variety of physical processes. In static interpretations, the preferred 
direction is a spatial direction, and DP could represent the percolation of fluid through 
porous rock with a certain fraction of the channels blocked (De'Bell and Essam 1983b), 
crack propagation (Kertesz and Viscek 1980) or electric current in a diluted diode net- 
work (Redner and Brown 1981). In d^/namica/ interpretations, the preferred direction 
is time, and DP is modelled by a stochastic cellular automaton (Kinzel 1985) in which 
all lattice sites evolve simultaneously and the main interpretation is as an epidemic 
without immunisation (Harris 1974, Liggett 1985). The behaviour of these models 
is generally controlled by a single parameter p, which could be the probability that a 
channel is open or the infection probability depending on one's favourite interpretation. 

When p is smaller than a critical value Pc, the fluid docs not percolate through the rock 
(the epidemic dies out). Let P{p) be the probability that the wetted region percolates 
inflnitely far from the source (the ultimate survival probability in epidemic language) 
then one expects: 

P{p) OC (p - Pc)^, p p+. (1) 

DP type transitions are also encountered in many other situations, perhaps most promi- 
nently in Reggeon field theory (Grassberger and Sundemeyer 1978, Cardy and Sugar 
1980), chemical reactions (Schlogl 1972, Grassberger and dc la Torre 1979), in nu- 
merous models for heterogeneous catalysis and surface reactions (Ziff et al. 1986, 
Kohler and ben-Avraham 1991, Zhuo et al. 1993, Jensen 1994), self-organized critical- 
ity (Obukhov 1990 and Paczuski et al. 1994) and even galactic evolution (Schulman 
and Seiden 1982). This short and far from complete list clearly demonstrates that 
directed percolation is a problem which emerges in a diverse set of physical problems 
and therefore deserves a great deal of attention. 

In this paper we discuss series expansions for the percolation probability on directed 
square and honeycomb lattices. The earliest series expansion for the square bond 
problem was the eight terms calculated by Blease (1977). A great improvement was 
due to Baxter and Guttmann (1988) who extended this series to 41 terms. For the 
honeycomb bond problem Onody (1990) obtained a 13 term series and for the square 
site problem the longest series of 16 terms is due to Onody and Neves (1992) improv- 
ing the previous record of 10 terms held by De'Bell and Essam (1983a). Using the 
finite-lattice method pioneered in this context by Baxter and Guttmann (1988) we 
have extended these series to 54 terms for the square bond problem, 37 terms for the 
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square site problem and 36 terms for the honeycomb bond problem. The percolation 
probability for the honeycomb site problem is related very simply to that of the square 
site problem, P"^{p) = pSQ{p^) (Dhar et al. 1982, Essam and De'Bell 1982). Note 
also that bond percolation on the honeycomb lattice may be viewed as site-bond per- 
colation on the square lattice (Essam and De'Bell 1982). In passing, we note that long 
series have been obtained for the moments of the pair connectedness for the site and 
bond problems on square and triangular lattices (Essam et al. 1986, 1988). 



2 The finite- lattice method 



We wish to calculate the series expansion of the percolation probability on square 
and honeycomb lattices oriented as in Figure |l|. We shall consider both site and 
bond percolation on these lattices. In site (bond) percolation each site (bond) is 
independently present with probability p and absent with probability q = 1 — p. Two 
sites are connected if one can find a path passing through occupied sites (bonds) only, 
while always following the allowed directions. For an infinite system, when q is less 
than a critical value qc, there is an infinite cluster spanning the lattice. The order 
parameter of the system is the percolation probability P{q), i.e., the probability that a 
given site belongs to the infinite cluster. Note that a path passing through a given site 
can only lead to the sites shown in Figure ^ below the origin O. This naturally leads 
one to consider a finite-lattice approximation to P{q), namely the probability PN^q) 
that the origin is connected to at least one site in the iV'th row. Pn{ci) is a polynomial 
in q with integer coefficients and a maximal order determined by the total number of 
sites (bonds) that may be present on the finite lattice. 




Figure 1: The directed square and honeycomb lattices with orientation given by the 
arrows. The rows are labelled as according to the text. 
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It has been proved (Bousquet-Melou 1995), for all the problems considered in this 
paper, that the polynomials Pn{<i) have a formal limit in the algebra of formal power 
series in the variable q, and therefore P{q) = lim^v^oo -P/v(q')- In all cases one finds 
that PnIq) converges to P{q) in such a way that the first (or — 1 depending on 
the specific problem) terms of the polynomials Pn{<i) coincide with those of P{q)- 

2.1 Specification of the models 

In order to calculate -Piv(g) we associate a state with each site, such that = 1 
if site j is connected to the A^'th row and (Xj = —1 otherwise. We shall often write 
+/— for simplicity. Note that a site can be in state —1 even though, in the case of site 
percolation, it is itself occupied, or, in the case of bond percolation, bonds emanating 
from the site are present. Let l,r denote the sites below t as in Figure |l]. We then 
define the triangle weight function W{at\<Ji,ar) as the probability that the top site t 
of the triangle is in state at, given that the lower sites I to the left and r to the right 
are in states cx; and cXr, respectively. One can then prove (Bidaux and Forgacs 1984, 
Baxter and Guttmann 1988) that 

PNiq) = T.I[Wi^tWi,ar), (2) 

{a} t 

where the product is over all sites j of the lattice above the A^'th row. The sum is over 
all values ±1 of each at, other than the topmost spin ai which always takes the value 
+1. The spins in the A^'th row are fixed to be +1. In short Pn{<i) is calculated as the 
sum over all possible configurations of the probability of each individual configuration. 

The weights W are listed in Table |. Obviously, W{—\ai,ar) = 1 — W{+\ai,ar). 
The remaining weights are easily calculated by considering the various possible ar- 
rangements of states and bonds. W{+\—,—) = because the top site is connected 
to the A^'th row if and only if at least one of the neighbours is connected. Let us 
next look at the remaining square bond weights. iy(-|-|+,+) = 1 — because the 
only bond configuration not allowed is both bonds absent which has probability g^. 
Finally, ) = ,+) = 1 — q because the bond connecting the two + 

states has to be present, which happens with probability p = 1 — q, and the other 
bond can be either present or absent. For the honeycomb bond problem we find that 
ai,ar) = (1 - q)W^'^{+\ai,ar) because if the top state is +1 the vertical 
bond has to be present. Note that one can think of the honeycomb bond problem as 
site-bond percolation on the directed square lattice where both sites and bonds are 
present with equal probability (Essam and De'Bell 1982). For the square site problem 
the weights are a little simpler since a site can be in state +1 only if it is present and 
W picks up only the probabihty of the top state, therefore — , — ) = as before 
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Problem 




W{+\+,-) = Wi+\-,+) 


W{+\-,-) 


SQ-bond 




l-q 





HC-bond 




{l-q? 





SQ-site 


l-q 


l-q 





HC-site 


(1 - 


{l-qf 






Table I: The triangle weight functions for the various directed percolation problems. 
Generally one has W{—\ai^ a^) = 1 — lT^(+|cr/, cr,.). 

and all the other weights with a +1 top state are equal to 1 — g. The honeycomb site 
weights are derived from the square site weights in the same manner as for the bond 
case. Note that it is costumary to assume in site percolation problems that the origin 
is present with probability 1. 

For the square and honeycomb site problems we therefore find: 

PN{q) = Y.X{^{^Muar) = Y.Wo{'^\(T2,(T^)\{W{at\auar) 

{<x} t {<x} t 

The weights W{at\cri,ar) are those of Table | and Wq is the weight of the top-most 
triangle. It is clear from Table | that for the site problem W^^{q) = W^^{2q — q^). 
Since the 'top' weights are 1 for the square site problem and 1 — q for the honeycomb 
site problem we find that 

^r(?) = (i-9)^;?'^(i-(i (3) 

which is essentially the relation mentioned in the Introduction, derived from the work 
of Dhar et al. (1982) by Essam and De'Bell (1982). 

2.2 Series expansion algorithm 

For small it is quite easy to calculate -P/v(q') by hand, but for larger one obviously 
has to resort to computer algorithms. The algorithms are basically implementations of 
a transfer matrix method. From Eq. @) we see that the evaluation of Pniq) involves 
only local 'interactions' since the weights involve only three neighbouring sites. The 
sum over all configurations can therefore be performed by moving a boundary line 
through the lattice. At any given stage this line cuts through a number of, say m, 
lattice sites thus leading to a total of 2"^ possible configurations along this line. Any 
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Figure 2: Part of the directed square lattice with the present boundary indicated by open 
circles. All weigths to the left of this boundary have been summed up. The weigth of the 
triangle given by (xr, x',., Xr~i) is picked up by moving the boundary from Xr to and 
updating the associated polynomials according to Eq. ^. 

configuration along the line is trivially represented as a binary number by letting the 
r'th bit of the number equal {ar + l)/2. For each configuration along the boundary 
line one maintains a (truncated) polynomial which equals the sum of the product of 
weights over all possible states on the side of the boundary already traversed. The 
boundary is moved through the lattice one site at a time. In Figure |^ we show how 
the boundary is moved in order to pick up the weight associated with a given triangle 
at position r along the boundary line. Let 5*0 = ( 

configuration of sites along the boundary with at position r and similarly SI = 
{xi, Xr-i, 1, Xr+1, Xm) the Configuration with 1 at position r. Then in moving the 
r'th site from the bottom left to the top of the triangle we see that the polynomials 
associated with these configurations are updated as 

P{SO) = W{0\0,Xr-i)P{SO) + W{0\l,Xr-l)P{Sl), 

(4) 

P{S1) = W{l\0,Xr^i)P{SO) + W{l\l,Xr-l)P{Sl). 

The calculation of Pn{(i) by this method is limited by memory, since one needs stor- 
age for 2^^^ boundary configurations. To alleviate this problem one can introduce 
a cut into the lattice, fix the states on this cut, evaluate the lattice sum P^{q) for 
each configuration C of the cut, and finally get Pn{(i) = He Pn^.^) the sum over 
all configurations of the cut. By placing the cut appropriately the growth in memory 
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requirements can be reduced to 2^/^. Obviously the finite- lattice calculation for dif- 
ferent configurations of the cut are independent of one another and these algorithms 
are therefore perfectly suited to take full advantage of modern massively parallel com- 
puters. In the following section we give a few more details of the algorithms we have 
used. 

2.2.1 The bond problem algorithm 

A very efficient algorithm was devised by Baxter for the square bond problem (Baxter 
and Guttmann 1988). All we had to do for the present work was basically to parallelize 
the algorithm in order to fully utilize the Intel Paragon at Melbourne University. The 
algorithm is based on an ingenious transformation of the square bond problem onto a 
honeycomb lattice. This is done by noting that the square bond weights can be written 
as 

W{at\cri,ar) = fi(^t,crm)gicrh(Tm)g{(^r,(Tm) (5) 

(7m=±l 

where 

/(+, +) = -!, /(+, -) = /(-, +) = 1, /(-, -) = 

9i+, +) = g, 9i+, -) = 9i-, +) = 9i-, -) = 1- (6) 

This means that if we replace each upwards pointing triangle in Figure |l| by a three- 
pointed star, arriving at the honeycomb lattice of Figure |, then Pn{(i) can be calcu- 
lated from this lattice by assigning weights /(cxj, aj) to vertical edges and g{cri, aj) 
to non-vertical edges. 

A cut of length L is introduced along the line RS in Figure |^ and the transfer matrix 
technique is used to build up the lattice to the left of RSO, starting from XR and 
working upvards to OS. The lattice is symmetrical around the central axis RSO and 
one can therefore obtain the lattice sum for the whole lattice by forming the sum of 
the squares of each boundary line polynomial. After this operation, the whole lattice 
is summed except for the edges on the center-line. So finally one has to multiply the 
(squared) boundary line polynomials by the weights of these edges. This is where 
the great advantage of the transformation becomes clear. Because /(— , — ) = we 
need never consider configurations of the cut (or parts of the boundary line in the 
vertical position) which have any (— , — ) edge. This basically means that the number 
of configurations of the cut which contribute to -P/v(q') grow only like 3^^^ rather than 
the usually expected 2^. The transformation thus provides us with an exponentially 
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Figure 3: The transformed lattice used in the square bond algorithnn. The sites nnarked 
with the full circles on the cut-line SR are fixed. 

faster algorithm. Likewise, as parts of the boundary hue enter the vertical position no 
(— , — ) edges need to be considered which leads to a significant reduction in the length 
of the cut for a given amount of memory. The memory requirement for the algorithm 
is governed by the maximal extent of the boundary line, which is at XR, and hence 
grows hke 2^/^~^ With this algorithm we calculated PN{q) for N < 39. Since the 
integer coefficients occuring in the series expansion become very large the calculation 
was performed using modular arithmetic (see, for example, Knuth 1969). Each run, 
using a different modulus, took approximately 24 hours using 50 nodes on an Intel 
Paragon. 

Virtually the same algorithm can be used for the honeycomb bond or site-bond square 
problem except that the / weights have to be replaced by 



/(+,+) = -(l-g), /(+,-) = /(- +) = l-g, /{-,-)= q. (7) 

Since /(—,—) no longer equals 0, obviously the great advange of the original trans- 
formation vanishes and the number of configurations of the cut grow hke 2^. For this 
reason we had to stop calculating Pn{q) ai N — 33, where each modulus required 
about 32 hours of CPU time using 50 nodes. 
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2.2.2 The site problem algorithm 

For the site problem the growth in memory can be limited to 2^^^^^ by introducing a 
cut across the lattice at row N/2. The upper part of the lattice is built up first by the 
transfer matrix technique, yielding a partial lattice sum and then the lower part 
is done. The total lattice sum for a given cut P^ is simply the product of these, 
i.e., P^ = P^P^j . Again Pn^q) is the sum over all configurations of the cut. It might 
seem that the number of cuts grow as 2^/^. Substantial simplifications can however 
be obtained. Note first of all that there is symmetry around the central vertical line 
which basically reduces the number of cut-configurations by a factor of 2. A more 
subtle means of reducing the number of cuts is obtained as follows: Since all triangle 
weights with at least one + on the bottom are the same, it follows that for any two 
configurations C and C which can be turned into one another by changing any number 

of +'s to — 's without adding or removing any ( ) sequences, P^ = P^' . This means, 

for example, that for any cut C without ( ) occurences, equals the partial sum 

of the all +'s cut. It is possible to use this property to perform the lower lattice sum 
simultaneously for many cuts. As an example consider cuts starting with ++ and — h 
but otherwise the same. The upper part is the same and the lower part is also the 
same except for the weight of the left-most triangle on the cut. By considering the 
various possibilities when moving the boundary line across this point, one can easily 
see that the two configurations can be summed simultaneously, i.e., the — h cut can 
be made as part of the -|--|- cut. We calculated Pn{(i) for < 32 which took about 48 
hours for each modulus with A^ = 32 using 50 nodes. 

We also calculated the series expansion for the honeycomb site problem up to P32(q'). 
Although we know the exact relation between the two site problems, Eq. (j^), this 
calculation provides us with an extra check of the algorithm and the extrapolation 
formulas we shall discuss presently. 

3 Extrapolation of the series 

As mentioned Pn{q) will generally agree with the series for P{q) up to some order de- 
termined by A^. For the square bond problem the coefficients of PNiq) = Sm>o iN,m<f^ 
agree with those of P{q) = Y.m>o C'mQ"^ to order A^. Baxter and Guttmann (1988) found 
that the series for P{q) can be extended considerably by determining the correction 
terms to Pn{(i)- Let us look at 

PN-PN+l = q''Y.(l''dN,r (8) 

r>l 
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then we shall call d^^r = 0'N,N+r — 0'N+i,N+r the r'th correction term. Obviously if one 
can find formulas for dN,r for allr < K then one can use the series coefficients of Pn{ci) 
to extend the series for P{q) to order N + K since 



k 

(iN+k — (J'N,N+k ~ X] ^N+k 
m=l 



■—m,m 



(9) 



for all k < K. That this method can be very efficient was clearly demonstrated by 
Baxter and Guttmann, who identified the first twelve correction terms, and used P29{(l) 
to extend the series for P{q) to 41 terms. To really appreciate this advance one should 
bear in mind that the time it takes to calculate Pn{(i) grows exponentially with N, 
so a direct calculation correct to the same order would have taken years rather than 
days. In the following we will give details of the correction terms for the various cases. 

3.1 The square bond case 

The first correction term for the square bond case is given by the Catalan numbers 



a result which was proved (Bousquet-Melou 1995) by noting that the correction term 
arise from compact bond animals of directed height N and perimeter N + 1. The 
second correction term 



was also calculated exactly recently (Bousquet-Melou 1995). As noted by Baxter and 
Guttmann (1988) the higher-order correction terms d^^r can be expressed as rational 
functions of the Catalan numbers. We have found that dN,r always can be written in 
the form 



where m — max(0,r — 4 — 2k). These formulas hold for all available N, provided 
that only Catalan numbers Cm with m > are involved. As noted by Baxter and 
Guttmann it is also true for m = —1 provided one 'defines' c_i — —1 (there was a 
misprint at this point in the original article) . Thus the extrapolation formulas are true 
for > r — 4. For r < 15 the coefficients Ar^k and Brj are either integers or fractions 



dN,i ^ cn ^ {2N\)/N\{N + ly. 



(10) 



dN,2 — '2cn — Cn+1 



(11) 




(12) 
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Table III: The coefficients a„ in the series expansion of P{q) = Y,n>o(^n(f' ^or directed 
bond percolation on the square lattice. 

with small (2 or 5) denominators. Note that there are various relations between the 
Catalan numbers so there are infinitely many ways of writing (JT^). For several of the 
correction formulas the general form adopted in this paper is slightly different from 
that of Baxter and Guttmann (1988) who tried whereever possible to choose a form 
involving only integers. The trade-off for having a general expression for the correction 
terms is that more rational fractions become involved. However with the proliferation 
of powerful mathematical packages such as MAPLE and MATHEMATICA this trade- 
off is well worth while. In Table |T| we have listed the coefficients and -B^j for 
r < 15. Using these extrapolarion formulas and the series for P^Q{q) we have extended 
the series for P{q) to the 54 terms given in Table ^TT[ 
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3.2 The square site case 



Inspired by the success of the extrapolation procedure for the square bond problem one 
might hope for similar success for other problems. And indeed one can find several of 
the correction terms for the square site problem, although the success is less spectacular 
as one is restricted to the first six correction terms. The first correction term dN,i was 
identified by Onody and Neves (1992) and since computed exactly by Bousquet-Melou 
(1995) 

d (3iV)! 

^^'^ - N\{2N + l)V ^^^^ 

This expression for the correction term was identified by Onody and Neves as the 
number of ways of inserting n — 4 sheets through a ball having n vertices on its surface 
such that pairs of sheets meet only on surface curves joining vertices! While this is true, 
a more useful and pertinent interpretation can be given. Viennot (1994) has pointed 
out that this is just the expression for the number of ternary trees of n vertices, which 
in turn is isomorphic to the number of diagonally convex directed animals (Delest and 
Fedou 1989). It is the identification between these animals and the first correction 
term that has been proved by Bousquet-Melou. She also proved our formula for the 
second correction term. 

As in the square bond case we can express higher correction terms as a function of 
dN,i- Again, there are infinitely many ways of expressing the formulas for the correction 
terms, one of which is 

dN = Y.CrA ■ dN,l + E i^Br,, + A,,)^^,! (M) 
i=2 \ ^ / j=l 

where K = N — r + j for r < 4 and K = N — r — 1+j for r > 5. These formulas are 
correct up to r = 6, whenever N > r. The coefficients are listed in Table ^V\. These 
formulas allowed us to extend the series for P{q) by an additional six terms to a total 
37 terms listed in Table 0. 

3.3 The honeycomb bond case 

For bond percolation on the directed honeycomb lattice Bousquet-Melou (1995) proved 
that the generating function / = J2N>idN,it^~^ of the first correction term dN,i is 
characterized by the algebraic equation 
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cr 



E 

.o 



r/j 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



3 

-12 



-31 
11 

-73| 
23| 
18i 



4" 

-353 
-199 
-123 
249 i 
904 
-396i 

6 



5 

-1050 
-5458| 
-2965 1 
-5302i 

000 2014 

—AqqIM. 

^■^■^256 

-8771 if 
152511 
392 



6 

-87304 
-71645 
-85165 
-7647|i 
-5605il§ 

134 11888. 

59049 
169 
1536 
107 



-1225 



39744^^2 
20245^ 
-8251^ 
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46 i 

-7i 

' 3 



-12^ 



-277 



12 
s 11 
^12 

7^ 



Br 



126 
120| 
150^ 
948^ 
862 i 
-149^ 




225 
1388 1 
1113 
267| 
-404 i§ 

-1 CO/l 353035 
"^"'^^ 559872 

-7324^ 



362 



192 

112 



15372 
13474 1 
23201 1 
2194^ 
3084 Ifi 

POP 17965 
-888 19683 
CCQO 6619423 
oo^a 30233088 

31048 1|§ 
10870 |I 
-2273l§ 




Cr 
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22 



1679616 
10969 



-12 
11 



15552 
385 



56 
-36 



90699264 
2734421 
5038848 
30595 



19 



46656 
20|||i 



'+-' 

o 

Q. 
(TJ 
+-' 



-o 



■M 
C 

cu 



O 
u 



> 



n 




n 







1 


19 


-92459524 

C L J C O ^ 


1 


Q 


20 


-298142956 


2 


-1 


21 


-922424269 


3 


-3 


22 


-3098690837 


4 


-8 


23 


-90429371 79 


5 


-21 


24 


-34187149573 


6 


-56 


25 


-79544646085 


7 


-154 


26 


-439149878359 


8 


-434 


27 


-313237196088 


9 


-1252 


28 


-7786443675714 


10 


-3675 


29 


16637473844344 


11 


-10954 


30 


-207593240544002 


12 


-33044 


31 


973714665769453 


13 


-100676 


32 


-7311741153076579 


14 


-309569 


33 


43345744201832502 


15 


-957424 


34 


-292472879532946388 


16 


-2987846 


35 


1867850225746155582 


17 


-9330274 


36 


-12389925641797917900 


18 


-29522921 


37 


81441868912809214904 



Table V: The coefficients a„ in the series expansion of P{q) = I]„>o'^n5'" for directed 
site percolation on the square lattice. 

f = l + tf + lii7 + t)f + f). (15) 
The higher-order correction terms are given by the formulas 

2r 

dN-2r+2+k,l + y^(A^i?-r.j + Arj)dN-2r+2+j,l (16) 

i=i 

which we find to be correct for r < 4 and > 2r — 1. The coefficients are listed in 
Table |V| apart from the Z?'s since the only non-zero ones are D41 = —157281/5 and 
-04,2 = 1744273/5. The final 36 terms series for P{q) is given in Table |VI1| . 



k=l 




+ Cr. 
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r/j 2 3 4 


2 3 4 


2 3 4 


1 _8A io954i -2773464^ 

2 -67 i§ 87495 |i -24330909 iff 

3 -3f 183llf 415040^ 

4 593^ -276663| 

5 66i§ -56787| 

6 7| -8377f 

7 -646 i 

8 -315f 


14i -5696|i 814389^^ 

54§ -48502^ 7805538i^ 

4-^ — 2Q2— 'igiO'l'i 
^20 ^^^20 *^^^^° 2000 

-616^ 98675 II 
23 i 31375 g 
5i 9476i 

-187^ 


-12i| 1753^ -86965^ 
-141^ 19295^ _1118660i§§ 
1213 -1663ii 224372^ 



Table VI: The coefficients A^j, Bj.j and Crj in the extrapolation formulas for the hon- 
eycomb bond problem. 



n 




n 


(In 





1 


19 


-1103369168956 


1 


-1 


20 


-5771541600014 


2 


-4 


21 


-31153472926184 


3 


-12 


22 


-160153702442390 


4 


-45 


23 


-907425183546587 


5 


-188 


24 


-4317291410619157 


6 


-835 


25 


-28433248376749141 


7 


-3849 


26 


-99125481158184567 


8 


-18242 


27 


-1076035285073833314 


9 


-88265 


28 


-238091850291444337 


10 


-434295 


29 


-58631611223043405378 


11 


-2165198 


30 


279283045229982597450 


12 


-10915089 


31 


-4730770444199592196256 


13 


-55534781 


32 


40182669640102878093220 


14 


-284708699 


33 


-480633574529182764438221 


15 


-1470350760 


34 


4852667371105928333619923 


16 


-7628363273 


35 


-53829647651783620888423836 


17 


-39878267745 


36 


574209696129704803372604206 


18 


-208458228964 


37 





Table VII: The coefficients a„ in the series expansion of P{q) — J2n>oC'n(l^ for directed 
bond percolation on the honeycomb lattice. 



15 



4 Analysis of the series 



We expect that the series for the percolation probabihty behaves hke 

P{q) ~ A{1 - q/q,Y[l + a^{l - q/q,)"" + ...], (17) 

where A is the critical amplitude, A the leading confluent exponent and the . . . repre- 
sents higher order correction terms. By universality we expect 13 to be the same for all 
the percolation problems studied in this paper and we will argue that the dominant 
correction term is analytic, i.e., A = 1. 

In the following sections we present the results of our analysis of the series which 
include accurate estimates for the critical parameters gc, A and A. For the most 
part the best results are obtained using Dlog Fade (or in some cases just ordinary 
Fade) approximants. A comprehensive review of these and other techniques for series 
analysis may be found in Guttmann (1989). 



4.1 qc and (5 

In Table |VlTT| we show the Dlog Fade approximants to the percolation probability series 
for bond percolation on the directed square lattice. The defective approximants, those 
for which there is a spurious singularity on the positive real axis closer to the origin than 
the physical critical point, are marked with an asterisk. The overwhelming majority 
of the approximants cluster around the values qc = 0.3552994 and P = 0.27643. As 
always in this type of analysis it is very difficult to accurately judge the true errors of 
the estimates of the critical parameters, however we venture to say that the critical 
parameters lie in the ranges, qc = 0.3552994(10) and /? = 0.27643(10), where the figures 
in parenthesis indicate the estimated error on the last digits. The other remarkable 



feature of Table |VIII| is that surprisingly many of the high-order approximants are 
defective. 



The results of the analysis of the series for the square site problem are listed in Table |IX. 
In this case there is a marked upward drift in the estimates for both qc and (3 and the 
estimates do not settle down to definite values. It does however seem likely that the 
true critical parameters lie within the estimates: qc = 0.294515(5) and (3 = 0.2763(3). 

The analysis of the series for the honeycomb bond problem yields the results in Table ^. 
Again we see an upward drift in the estimates for both qc and (3 though the estimates 
are somewhat more stable than in the previous case. It seems likely that the true 
critical parameters lie within the estimates: qc = 0.177143(2) and (3 = 0.2763(2). 

Finally we analysed the series for the honeycomb site problem, with the results tabu- 
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N 


[N-1,N] 


[N,N] 


[N+1,N] 






Qc 




p 




Qc 




/? 




Qc 






11 





3553000 





27645 





3553030 





27653 





3553023 





27651 


12 





3553016 





27649 





3553011 





27648 





3552997 





27644 


13 





3553028* 





27652* 





3553004 





27646 





3553000 





27645 


14 





3552994 





27643 





3552972 





27634 





3552995 





27643 


15 





3552991 





27642 





3552994 





27643 





3552994 





27643 


16 





3552994 





27643 





3552994 





27643 





3552994 





27643 


17 





3552994 





27643 





3552994 





27643 





3552997* 





27644* 


18 





3552994 





27643 





3552992 





27642 





3552983 





27632 


19 





3553002* 





27643* 





3552991 





27641 





3552996* 





27644* 


20 





3552994 





27643 





3552994 





27643 





3552994 





27643 


21 





3552994 





27643 





3552994 





27643 





3552994* 





27643* 


22 





3552994* 





27643* 





3552994* 





27643* 





3552994* 





27643* 


23 





3552994* 





27643* 





3552994 





27643 





3552994 





27643 


24 





3552993* 





27643* 





3552993* 





27643* 





3552995 





27644 


25 





3552993* 





27643* 





3552997 





27645 





3552995* 





27644* 


26 





3552991* 





27643* 





3552990* 





27643* 





3552986* 





27647* 


27 





3552993* 





27643* 



















Table VIII: Dlog Pade approximants to the percolation series for directed bond percolation 
on the square lattice. 



N 


[N-1,N] 


[N,N] 




[N+ 


1,N] 








Qc 








Qc 




/3 




Qc 






5 





2939337 





26881 





2943291 





27266 





2942979 





27228 


6 





2942670 





27190 





2943175 





27252 





2942699* 





27199* 


7 





2944168 





27393 





2944521 





27453 





2944777 





27502 


8 





2945135 





27585 





2944742 





27495 





2944794* 





27505* 


9 





2944599 





27465 





2944720 





27490 





2944739 





27494 


10 





2944753 





27498 





2944656* 





27478* 





2944942 





27546 


11 





2945228 





27655 





2945156 





27623 





2945020 





27571 


12 





2945246* 





27662* 





2945060 





27586 





2945058 





27585 


13 





2945058 





27585 





2945061* 





27586* 





2945047 





27581 


14 





2945051 





27578 





2945051 





27582 


* 




* 




15 





2945056 





27584 





2945047* 





27581* 





2945032* 





27576* 


16 





2945069 





27589 





2945096 





27602 





2945089 





27598 


17 





2945090 





27599 





2945095 





27601 





2945113 





27612 


18 





2945134 





27625 





2945111 





27611 











Table IX: Dlog Pade approximants to the percolation series for directed site percolation 
on the square lattice. -ly 



N 


[N-1,N] 


[N,N] 


[N+1,N] 






Qc 




/? 

r 




Qc 




1 




Qc 




p 


5 





1770229 





27331 





1770722 





27420 





1771131 





27507 


6 





1771195 





27523 





1770967 





27469 





1771067 





27493 


7 





1771087 





27498 





1771161 





27517 





1771270 





27552 


8 





1771320 





27572 





1770209* 





27662* 





1771414 





27612 


9 





1771480 





27647 





1771294 





27559 





1771369 





27591 


10 





1771391 





27601 





1771352 





27584 





1771356 





27585 


11 





1771357 





27586 





1771344* 





27580* 





1771399 





27609 


12 





1771412 





27619 





1771381 





27598 





1771395 





27606 


13 





1771402 





27612 





1771411 





27618 





1771403 





27612 


14 





1771406 





27614 





mum 





27641 





1771403* 





27612* 


15 





1771405 





27613 





1771408 





27616 





1771429 





27636 


16 





1771390* 





27605* 





1771415 





27622 





1771419 





27625 


17 





1771422 





27629 





1771418 





27624 





1771418 





27624 


18 





1771418 





27624 



















Table X: Dlog Pade approximants to the percolation series for directed bond percolation 
on the honeycomb lattice. 



N 


[N-1,N] 


[N,N] 


[N+1,N] 










p 




Qc 




p 




Qc 




p 


5 





1598159 





27017 





1599573 





27265 





1599491 





27249 


6 





1599269 





27203 





1599516 





27254 





1599487* 





27248* 


7 





1600181 





27416 





1600409 





27485 





1600545 





27532 


8 





1600656 





27577 





1600476 





27507 





1599682 





27195 


9 





1600378 





27473 





1600457 





27501 





1600452 





27499 


10 





1600453 





27499 





1600456 





27501 





1600555 





27543 


11 





1600280* 





27462* 





1600711 





27640 





1600597 





27565 


12 





1600498 





27515 





1600630 





27585 





1600630 





27585 


13 





1600630 





27585 





1600630 





27585 





1600622 





27580 


14 





1600620 





27579 





1600625 





27582 





1600636 





27589 


15 





1600630 





27585 





1600622* 





27580* 





1600391* 





27665* 


16 





1600641 





27593 





1600656 





27606 





1600647 





27597 


17 





1600650 





27600 





1600655 





27604 





1600662 





27611 


18 





1600688 





27642 





1600662 





27611 











Table XI: Dlog Pade approximants to the percolation series for directed site percolation 
on the honeycomb lattice. 
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lated in Table 0. As in the square site case thereis a very pronounced upward drift 
in the estimates for both and 13. It seems hkely that the true critical parameters lie 
within the estimates: qc = 0.160067(5) and /? = 0.2763(4). We note that the expected 
relation between the values of qc for the square site and honeycomb site problems, 
gSQ _ 2q^^ — {q^^y, clearly is fulfilled by the estimates. This inspires some con- 
fidence in the appropriateness of our extrapolation method in general and our error 
estimates in particular. 



4.2 The critical amplitudes 

From the leading critical behaviour, P{q) ~ A{1 — q/qc)^, it follows that {qc — 
q)P^^^^\q=g^ ~ A'^/^qc. So by forming the series for G{q) = {qc — q)P~^^^ we can 
estimate the critical amplitude A from Fade approximants to G evaluated at qc. The 
prodecure works well but requires knowledge of both qc and (3. For the square bond 
series we know both qc and /3 very accurately, and we estimated A using values of qc 
between 0.355299 and 0.3553 and values of f5 ranging from 0.2764 to 0.2765. For each 
{qc, (3) pair we calculate A as the average over all [A^ + K, N] Fade approximants with 
A' = 0, ±1 and 2A^ + K > 45. The spread among the approximants is minimal for 
qc = 0.3552994, p = 0.27643 where A = 1.3291475(2). Allowing for values of qc and P 
within the full range we get A = 1.3292(5). 

For the square site series we used values of qc from 0.29451 to 0.29452 and (3 from 
0.2761 to 0.2765 averaging over Fade approximants with 2A^ + K > 27. In this case 
the spread is minimal for qc = 0.294515, (3 = 0.2763 with A = 1.425164(5). Again 
allowing for a wider choice of critical parameters we estimate that A = 1.425(1). 

For the honeycomb bond series we restricted qc to lie between 0.177138 and 0.177148 
and (3 between 0.2761 to 0.2765 using all approximants with 2A^ + K > 26. The 
minimal spread occurs at qc = 0.177143, /3 = 0.27635 where A = 1.10607(2). A wider 
choice for qc and f3 leads to the estimate A = 1.106(1). 

Finally in the honeycomb site case we used values of qc in the range 0.160065 to 
0.160075 and P from 0.2761 to 0.2765 using all approximants with 2N + K > 27. 
The minimal spread occurs when = 0.160069, /? = 0.2764 where A = 1.16779(2). 
With the wider choice of critical parameters we estimate that A = 1.167(1). The 
exact relation, Eq. (|^), between the square and honeycomb site problems means that 
there is a simple relation between the amplitudes in the two cases. First note that, 
^""(l - q/QcH)" ~ P^(g) = (1 - g)P^(2g - g^) ^ (l _ g)(l _ (2g - q^)/qc,sy. Since, 
QcS = '2qc,H - qIh, we find that, (1 - (2g - q^)/qc,sY = [{qc,H - g)(2 - qc,H - q)/qc,sY, 
and therefore: = (1 - q^uW ' 2qc,H)qc,H I q^sf A^ = (1 - ge,H)(l - qlu/qcsfA'. 
Insertions of the various critical parameters shows that this relation is indeed satisfied 
by our amplitude estimates. 
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L 


SQ bond 


SQ site 


HC bond 


HC site 


1 


1.29661 


1.41614 


1 


11520 


1 


16740 


2 


1.31234 


1.39775 


1 


12002 


1 


16579 


3 


1.31114 


1.39989 


1 


12001 


1 


16607 


4 


1.31218 


1.37739 


1 


11952 


1 


16564 


5 


1.31098 


1.39359 


1 


11750 


1 


16546 


6 


1.31006 


1.39001 


1 


11808 


1 


16521 


7 


1.32566 


1.39889 


1 


11856 


1 


16486 


8 


1.30916 


1.39582 


1 


11929 


1 


16537 


9 


1.31322 


1.39162 


1 


11929 


1 


16534 


10 


1.31122 


1.39449 


1 


11780 


1 


16508 


11 


1.31195 


1.40570 


1 


12056 


1 


16578 


12 


1.31228 


1.40306 


1 


12435 


1 


16462 



Table XII: Critical amplitudes A for the four percolation problems obtained by using 
the method of Liu and Fisher. The estimates were calculated by averaging over various 
inhomogeneous differential approximants of order L. 

A second method, proposed by Liu and Fisher (1989), for calculating critical ampli- 
tudes starts by assuming the functional form P{q) ~ A{q){l — q/qc)^ + B{q). One 
then transforms this function into g{q) = (1 — q/qc)^^P{q) ~ A{q) + B{q){l — q/qc)'^- 
The required amplitude is now the background term in g{q), which can be obtained 



from inhomogeneous differential approximants (Guttmann (1989) p89). In Table pCII 
we have listed the estimates obtained by averaging over various first order differen- 
tial approximants using at least 40 terms of the series for the square bond case and 
at least 25 terms in the other cases. The critical parameters qc and f3, used in the 
transformation of the series, were the central values of the estimates from the previous 
section. This method generally yields slightly lower estimates for the amplitudes and 
the spread among the approximants is much larger than in the first method. 



4.3 The confluent exponent 

We studied the series using two different methods in order to estimate the value of 
the confluent exponent. In the first method, due to Baker and Hunter (1973), one 
transforms the function, P{q) = J27=i^ii^ ~ q/Qc)^^' = J2'^=o O'nq"' , into an auxiliary 
function with simple poles at 1/Aj. We first make the change of variable q = qc{l — e~^) 
and find, after multiplying the coefficient of by k\, the auxiliary function 
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N oo ^ A 

HO = E E ^^(^^0^= = E T3T7' (18) 

i=l k=0 i=l *^ 

which has poles at ^ = 1/Aj with residue —Ai/Xi. The great advantage of this method 
(when it works) is that one obtains simuhaneous estimates for many critical parame- 
ters, namely, /3, A, and the critical amplitude, while there is only one parameter, in 
the transformation. In Figure |^ we have plotted, respectively, f3 and A as a function 
of the transformation parameter Qc for various [N ± K, N] Fade approximants, with 
N > 25. For each 'guess' for qc we performed the Baker- Hunter transformation and 
located the numerically largest and next-largest poles, which are the estimates for the 
reciprocals of —(3 and — (/5 + A), respectively. The majority of the approximants have 
a very narrow crossing region close to qc = 0.3552996(3), with (3 = 0.27645(3) and 
A = 1.000(5). In Table pCIII| we have listed the estimates for f3, A and the corre- 
sponding critical amplitudes obtained from the Baker-Hunter transformed series with 
qc = 0.3552996. The results strongly suggest that the leading correction to scaling 
term is analytic. Furthermore we note that the estimates for the critical amplitudes 
fully agree with those obtained from the first method used in the previous section. 



N 


M 


f3 


A 


A 


A X aj\ 


22 


23 


0.27645 


1.32925 


1.00097 


1.03202 


23 


23 


0.27646 


1.32930 


1.00013 


1.03029 


24 


23 


0.27863 


1.32369 


0.98439 


1.01224 


23 


24 


0.27645 


1.32925 


1.00090 


1.03181 


24 


24 


0.27647 


1.32931 


0.99994 


1.02993 


25 


24 


0.27549 


1.33100 


1.01375 


1.05322 


24 


25 


0.27645 


1.32926 


1.00078 


1.03149 


25 


25 


0.27648 


1.32935 


0.99922 


1.02857 


26 


25 


0.27589 


1.33038 


1.00698 


1.04048 


25 


26 


0.27645 


1.32926 


1.00064 


1.03114 


26 


26 


0.27649 


1.32936 


0.99906 


1.02826 


27 


26 


0.27617 


1.32992 


1.00305 


1.03410 


26 


27 


0.27645 


1.32928 


1.00037 


1.03052 


27 


27 


0.27649 


1.32936 


0.99911 


1.02836 



Table XIII: Estimates for the critical exponent (3, critical amplitude A, confluent exponent 
A, and confluent amplitude AxaA, obtained from [A^, M] Pade approximants to the Baker- 
Hunter transformed square bond series with qc = 0.3552996. 
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Figure 4: The critical exponent j3 and confluent exponent A as a function of the parameter 
qc in the Baker-Hunter transformation for the square bond series. 

In the second method, due to Adler et al. (1981), one studies Dlog Fade approximants 
to the function G{q) = PP{q) + {qc — (l)(iP{(l)/dq. The logarithmic derivative to G{q) 
has a pole at qc with residue /3 + A. We evaluate the Dlog Fade approximants for a 
range of guesses for qc and p. For each such guess we thus find an estimate for A; for 
the correct value of qc and (3 one would expect to see a convergence region in {qc, (3, A)- 
space. In practice we always froze either or (3 and examined A as a function of the 
other parameter. Figure ^ shows, respectively, A as a function of qc with (3 = 0.27643 
and A as a function of (3 with = 0.3552994. This analysis clearly support A ^ 1, 
and thus that there is no sign of any non-analytic corrections to scaling. 

For the square site series the results from the Baker-Hunter transformation is less 
convincing, as there is no value of qc at which the various aprroximants cross. If we look 
closely at the approximants evaluated at qc = 0.294515 we find, generally speaking, 
that only the [A^ — 1, A^] approximants yield estimates of (3 close to the expected value 
with corresponding estimates for A consistent with an analytic correction. The method 
of Adler et al. confirms that A ~ 1. 

In the honeycomb bond case several of the approximants to the Baker-Hunter trans- 
formed series has a crossing for qc = 0.177144(1), P = 0.2767(1) and A = 0.89(2), 
though it should be noted that the scatter is quite large. When we analyse the series 
using the second method we find that, for qc and (3 close to the central values from the 
Dlog Fade analysis, a value of 1 for A is fully compatible with the results. 
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Figure 5: The confluent exponent A as a function of, respectively, the parameter qc (with 
p = 0.27643) and the parameter p (with = 0.3552994) using the method of Adier et 
al. (1981). 



5 Conclusion 



In this paper we have presented extended series for the percolation probability for 
site and bond percolation on the square and honeycomb lattices. The analysis of the 
series leads to improved estimates for the percolation threshold (particularly for the 
honeycomb bond problem) and the order parameter exponent (5. To summarise we 
estimate that 



= 0.3552994(10), 


P = 


-- 0.27643(10) 




1.3292(5) square bond problem. 


qc = 0.294515(5), 




= 0.2763(3) 


A = 


1.425(1) square site problem. 


qc = 0.177143(2), 




= 0.2763(2) 


A = 


1.106(1) honeycomb bond problem 


qc = 0.160067(5), 




= 0.2763(4) 


A^ 


1.167(1) honeycomb site problem. 



The estimates for q^ = 1 — Pc for the square bond and site problem are in excellent 
agreement with those obtained by Essam et al. (1986, 1988), qc = 0.355303(6) and 
qc = 0.29451(1), respectively. The estimates for j3 clearly show, as one would expect, 
that all the models studied in this paper belong to the same universality class. The 
value of P does not suggest any simple fraction. Indeed, around the central value for 
P (square bond) we find only four fractions with denominators less than 1500. They 
are: ^ = 0.276422 ■ ■ ^ = 0.276429 ■ ■ fi = 0.276433 ... and if = 0.276435 . . .. 
None of these are remotely compelling, and leave open the question as to why this 
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apparently simple problem has such an ugly exponent. This does seem to be a frequent 
characteristic of directed problems, as evidenced by the recent study of the longitudinal 
size exponent of square lattice directed animals (Conway and Guttmann 1994) in which 
it was found that i/y = 0.81722(5), a result which suggests no simple rational fraction. 
Finally we note that none of the series show any evidence of non-analytic confluent 
correction terms. This provides a hint that the model might be exactly solvable. 

Acknowledgements 

Financial support from the Australian Research Council is gratefully acknowledged. 

References 

Adler J, Moshe M and Privman V 1981 J. Phys. A:Math. Gen. 17 2233 

Baker G A and Hunter D L 1973 Phys. Rev. B 7 3377 

Baxter R J and Guttmann A J 1988 J. Phys. A:Math. Gen. 21 3193 

Bidaux R and Forgacs G 1984 J. Phys. A:Math. Gen. 17 1853 

Blease J 1977 J. Phys. CSolid State Phys. 10 917 

Bousquet-Melou M 1995 Percolation models and animals Preprint 

Broadbent S R and Hammersley J M 1957 Proc. Camp. Phil. Soc. 53 629 

Cardy J L and Sugar R L 1980 J. Phys. A:Math. Gen. 13 L423 

Conway A R and Guttmann A G 1994 J. Phys. A: Math. Gen. 27 7007 

De'Bell K and Essam J W 1983a J. Phys. A:Math. Gen. 16 3145 

De'BeU K and Essam J W 1983b J. Phys. A:Math. Gen. 16 3553 

Delest, M.-P and Fedou, J.M 1989 Exact formulas for fully diagonal compact animals 
Internal Report LaBRI 89-06, Univ. Bordeaux I, France 

Dhar D, Phani M K and Barma M 1982 J. Phys. A:Math. Gen. 15 L279 

Essam J W and De'BeU K 1982 J. Phys. A:Math. Gen. 15 L601 

Essam J W, De'BeU K, Adler J and Bhatti F M 1986 Phys. Rev. B 33 1982 

Essam J W, Guttmann A J and De'BeU K 1988 J. Phys. A:Math. Gen. 21 3815 

Grassberger P and Sundermeyer K 1978 Phys. Lett. 77B 220 



24 



Grassberger P and de la Torre A 1979 Ann. Phys. NY 122 373 

Guttmann A J 1989 Asymptotic Analysis of Power-Series Expansions Phase Transi- 
tions and Critical Phenomena vol 13, eds. C. Domb and J. L. Lebowitz (Academic 
Press, New York) . 

Harris T E 1974 Ann. Prob. 2 969 

Jensen I 1994 J. Phys. A:Math. Gen. 27 L61 

Kertesz J and Viscek T 1980 J. Phys. CSolid State Phys. 13 L343 
Kinzel W 1985 Z. Phys. B 58 229 

Knuth D E 1969 Seminumerical Algorithms (The Art of Computer Programming 2) 
(Addison- Wesley, Reading MA) 

Kohler J and bcn-Avraham D 1991 J. Phys. A:Math. Cen. 24 L621 
Liggett T M 1985 Interacting Particle Systems (Springer, New York) 
Liu A J and Fisher M E 1989 Physica 156A 35 
Obukhov S P 1990 Phys. Rev. Lett. 65 1395 
Onody R N 1990 J. Phys. AiMath. Gen. 23 L335 
Onody R N and Neves UPC 1992 J. Phys. AiMath. Gen. 25 6609 
Paczuski M, Maslov S and Bak P 1994 Europhys. Lett. 27 97 
Redner S and Brown A C 1981 J. Phys. A: Math. Gen. 14 L285 
Schlogl F 1972 Z. Phys. 252 147 

Schulman L S and Seiden P E 1982 J. Stat. Phys. 27 83 
Viennot X G 1994 Private communication 

Zhuo J, Redner S and Park H 1993 J. Phys. A:Math. Gen. 26 4197 
Ziff R M, Gulari E and Barshad Y 1986 Phys. Rev. Lett. 56 2553 



25 



r 





j 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 




1 


-2 


2 


2 


6 


28 


168 


1188 


9438 


81796 


758472 


7423832 


75951512 




2 


2 


-5 


13 


14 


40 


172 


940 


6022 


43224 


337852 


2821008 


24820084 




3 




5 


-29i 

Z 


58 


76 


230 


916 


4358 


22626 


116670 


490220 


-39210 




4 




-2 


-2 


-73 


218 


-23 


16720 


-238382 


5439166 


-54821200| 


1211505650| 


-8785667878 1 




5 






12 


24 


-360 


1185 


-2285 


66470 


-1329734 


13844727 1 


-3022267671 


2199183639| 




6 






-2 


-6 


174 


-1619 


3737 


20838 


-289156 


3246796 


-69778186 


503746552 




7 











-128 


889 


-7985 


23543 


-67332 


1 144926 i 


-20852206 1 


165626416| 


in ■ 


8 











82 


-242 


4320 


-35091 


56180 


355092 


-4943452 


44277656 


a^ ■ 

<^ : 


9 
10 










-15 



-15 
18 


-1334 
-104 


23956 
-8722 


-177605 
124506 


406391 1 
- 793053 1 


-1 134497 1 
1126018| 


12747113| 
3369765 i 




11 












-1 


243 


2820 


-56441 


675899 


-3956056 


7537547 


20 Sep 


12 















-46 


-1204 


18304 


-329150 


3495696 


-18307058 


13 














-2 


367 


-5702 


115282 


- 1913164 


18478108 


14 

















-44 


1636 


-25924 


715656 


-10857182 


15 
















2 


-250 


954 


-200132 


4463999 


16 



















8 


1036 


42000 


- 1380082 




17 


















-2 


-125 


-6857 


363300 


> ■ 


18 





















-12 


1064 


-95182 


^9509121 


19 




















2 


-163 


21166 


20 
21 
22 
23 
24 























12 
-2 



-2554 
116 
-12 
2 



■ 


25 

26 


























a: 


m 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


^ ■ 


1 


-2 


3 


^2 


13 


-90 


-745 


31152 


-485988 


10837718 


- 109826658 1 


2422283961 1 


-17572113989| 


iv:con 


2 
3 
4 
5 






4 


-2 


15 
-8 


-5 
4 


51 
-24 
16 


3566 
-19 



-57572 
24 
40 
-32 


854505 1 
-9414| 
81| 
-16 


-13937224| 
1011551 
-136| 
-72 


130944690 1 
-1480387| 
21604| 
-232 


arX. 


6 
7 






















64 


64 



Br 



lA. 

1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 



-6 
-12 



-31 
11 
-73| 
23| 
18| 



-353 
-199 
-123 
249 i 
904 
-396 i 




-1050 
-5458 i 
-2965 I 
-530 If 
293ii| 

193 

256 

183 



-499 



-8771256 
1525 f 
392 



-87304 
-71645 
-85165 

-7647|i 
-5605ii§ 

1 o /I 11888 
^'•'^ 59049 
-lOOr 169 
^^^O 1536 

39744 i2I 
20245^ 
-8251^ 




46i 
-7i 

' 3 



18 
-12| 
-277^ 
-812 



126 
120 1 
150^ 
948^ 
862 i 
-149^ 




-1524 



225 
1388 1 
1113 
267| 

-404 Ml 

353035 
559872 

-7324Tj^ 

2:H) 
"^"^ 192 



112 



15372 
13474 I 
23201 I 

2194A 

3084 fi 

J 17965 



' 19683 

gcno 6619423 
30233088 

31048 ifi 
10870§ 
-2273ii 




^36 



